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Abstract 

A new Majoron model is presented within the framework of the seesaw mechanism. 
Its Higgs sector consists of only doublet representations and the lepton-number viola- 
tion takes place at the same scale of the electroweak symmetry breaking. This model 
is different from the singlet- or triplet-Majoron model in several respects: it is free 
from the p-parameter constraint and it provides moderately fast neutrino decay, but 
the constraint from the stellar cooling of red giants is satisfied only with an imposed 
approximate symmetry. A r neutrino as heavy as 10 MeV is possible in this model 
despite various cosmological and astrophysical constraints. 



1. Introduction 

We present in this paper a new Majoron model. Its Higgs sector consists only of doublet 
representations under the standard electroweak gauge group G = SU(2)xU(l)y and the 
lepton number symmetry U(l)i breaks down spontaneously together with the breaking of 
G. The Majoron, the Nambu-Goldstone (NG) boson associated with the U(1)l breaking, 
is induced as a linear combination of only doublet Higgs bosons which is of course not the 
case in the singlet Majoron (SM) model or the triplet Majoron (TM) model [Q. In this 
respect, we refer to the present model as the doublet Majoron (DM) model. There is also 
a model in which the Majoron belongs to the doublet representation ||, but that version 
is now ruled out experimentally by the nonobservation of Z° — > invisible light scalars. In 
the present model the light neutrino masses come about through the seesaw mechanism [|J, 
while in the model of Ref. H] they are induced radiatively. 

One motivation of Majoron models in general is that they may induce neutrino decay 
through Majoron emission This decay can provide a way to evade the cosmological 
constraint on the neutrino mass m„, that is, m u < 100 eV or m v > a few GeV for a stable 
neutrino. In fact, this decay in both the TM and SM models is strongly suppressed and 
it turns out to be almost cosmologically irrelevant: in the TM model the decay only takes 
place through a loop diagram ||, while in the SM model it is suppressed at tree level by 
the fourth power of the seesaw factor {m v /M) (M is the mass scale for the right-handed 
gauge-singlet neutrino) ||. In the DM model, on the other hand, the decay occurs by the 
same mechanism as in the SM model, but the suppression of the seesaw factor is weaker; the 
power reduces to three instead of four. This means that neutrino decay in the DM model 
could well be cosmologically relevant. 

Some time ago, the possibility of a r neutrino of about 10 MeV in mass was proposed in 
a scenario of baryogenesis at the electroweak phase transition in the SM model . However 
it appears difficult to have such a mass range which also satisfies the cosmological constraint 
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unless one chooses an unnaturally small M. We aim to investigate the decay of such a 
neutrino more closely for a similar mass range using the DM model. We will also take 
into account the constraint of Majoron emission from the supernova 1987A ]§[] and that of 
nucleosynthesis P, |T0]. 



There is another motivation for looking into Majoron models. Although the seesaw 
mechanism is an attractive idea to account for the mass hierarchy of the neutrinos compared 
to the charged leptons J|, once it is combined with the observation that the anomalous 
baryon and lepton number B + L violating process is in thermal equilibrium at an early stage 



of the Universe |TTJ, it leads to a very strong bound on the neutrino masses |L2], [L3|]. If both 
the B + L violation from the anomalous process and the L violation from the gauge-singlet 
neutrino mass term are active together at a certain time during the Universe's evolution, 
any primordial excess in B or L is wiped out. To avoid this situation, a neutrino mass must 
be less than 1 eV as long as the anomalous process is fast enough up to a temperature of 
about 10 12 GeV This makes any neutrino an unsuitable candidate for a component of 
dark matter. 

One way to evade the above constraint is to make U(1)l an exact symmetry above the 
electroweak phase transition. In the DM model, the L symmetry is exact before spontaneous 
symmetry breaking even though the gauge-singlet neutral fermion has mass, and it breaks 
down at the electroweak scale, which is a natural consequence of its Higgs representation. 
Provided that this phase transition is first-order and that the anomalous process is suppressed 



after it [15], the L and B + L violating processes cannot coexist above or below the phase 



transition. 

2. The model 

We now describe the DM model. A basic observation that characterizes the DM model is 
that a Majorana mass for the gauge-singlet neutral fermions N a (a = 1,2,3 denotes genera- 
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tions) does not immediately imply L violation. Since N a belong to a different representation 
from the doublet leptons L a , they can have a lepton number different from that of L a . We as- 
sign zero lepton number to N a so that their Majorana mass terms do not violate L. Instead, 
we introduce a Higgs doublet Hi, different from the ordinary one H , and assign L — — 1 to 
it. (The hypercharge Y is the same for all the doublet Higgs bosons, Y = 1.) The Yukawa 
coupling between L a and N a is provided by Hi. The spontaneous breakdown of G, as H 
and Hi get vacuum expectation values, is accompanied by the violation of U(1)l. 

Actually, this simple extension of the minimal standard model to a Majoron model turns 
out not to be enough. The Majoron obtained this way has components in both H and Hi. 
As a result, the Majoron couples to the charged leptons through the hypercharge current jj 
as well as through the lepton current jjf. This coupling allows Majoron emission in Compton 
scattering and significantly contributes to the stellar cooling of red giants |5|, [16|]. To avoid 



this coupling, we introduce another Higgs doublet H2 and assign L = 1 to it. As we will 
see below, we can suppress the Majoron coupling to the charged leptons by a cancellation 
between Hi and H^. 

The Lagrangian is the same as in the standard model except for the part including the 
extra Higgs bosons and gauge-singlet neutral fermions. The corresponding part is 

C = J2\ D M 2 -v(H a ) 

a 

+ [iLl(D^)L a + iCl(D^)C a + iNl(d^)N a 

a 

+ E 



a,b 



HlT° ab (L T al a 2 C b ) + HlT l ab {Llia 2 N b ) + ^(A^ 2 ^) + (h.c.) 



(1) 



where all the lepton fields have been written in left-handed two-component notation; L a 
denotes the doublet leptons 

C a the singlet charged leptons e c , v c , and r c ; iV a the singlet neutral fermions; H = ir 2 H*; 



is the covariant derivative and a 11 = (1, — a); Note that Hi cannot have a Yukawa coupling 
because of its quantum number Y = 1 and L = 1. 

The Higgs potential which is invariant under G and U(l)^ consists of thirteen terms, i.e., 
three of the type (H^H a ), three of the type (H^H a ) 2 , three of the type (H^H a ) (H^H p) , three 
of the type (W a Hp)(H^H a ), and (HqHi)(HqH 2 ) plus its hermitian conjugate. For the last 
term, we have used the rephasing of H a to absorb a possible complex phase factor. Thus V 
has thirteen real parameters and has no explicit CP violation. For a wide range of parameters, 
V breaks GxU(1)l down to U(l) cm spontaneously with nontrivial real expectation values v a 
of H a . Let us write the scalar fields as 
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These thirteen parameters are phenomenologically constrained for the following reason. 
In the imaginary parts of the neutral sector of the Higgs bosons, there are two NG bosons 
associated with the violation of U(l)y and U(1)l. An appropriate linear combination is the 
one absorbed into Z°, and the other orthogonal combination remains massless and becomes 
the Majoron. The specific form of the linear combination for the Majoron in terms of ip a is 
obtained by working with the currents j£ and j'J, 

j L » = zJ2L(H a ){HlD»H a - (D*H a )*H a } + £ L(f)(p^f), (4) 

/ 

where L(H a ) and L(f) denote the lepton number assignment for H a and the fermion species 
/; is given by the same equation with Y(H a ) and Y(f), the hypercharges, instead of 
L(H a ) and L(f). We write the current-conservation equations using Eq. (|3]). They are 

d^Y,L{H a )v aVa = - 9 fe L{H a )vl\ d»Z» + -UmJ^OO + - (5) 

V2 cos 0iy \ a J V2 

0„0"5>«¥>« = - 9 fevA d»Z» + -^d,j Y "(f) + .... (6) 

a V2cOS#w V a J v2 

where 9yy is the Weinberg angle and (^) denotes the fermion component in each current. 
These equations describe the motion of the NG bosons; the right-hand sides represent the 



interactions with the other fields. The Majoron is the combination that does not have the 
linear term in its equation. We can thus simply use Eqs. fl5|) and (|6|) to extract out the 
equation for the Majoron. This procedure reveals that the Majoron has a coupling of the 
form g^eed^j Y ^ with the strength proportional to 

(7) 

a 

This coupling is highly constrained to be less than about 1CT 13 divided by the electron mass 
from the astrophysical consideration already mentioned |5|, 16 . 



The reason we employed H 2 is to cancel the Hi contribution in (0) [17]. Note that it 
vanishes for L(H\) = —L(H 2 ) and v\ = v 2 . Hence we parametrize V as 

V = J2\ a (HiH a -vl) 2 +J2Vap(HlH a -vl)(HlH p -v}) 

a a </3 

+ E Ca/3 [(HlH a )(HlH p ) - (HlH p )(HlH a ) 

a</3 

+£ [(HlH^HlHr) + {HlH Q ){H\H 2 ) - (H^H^Hi) - {H\H q ){hIh q )\ , (8) 

and further require it to be symmetric under the exchange Hi <-» H 2 . The parameters are 
then constrained: V\ = v 2 = vl, Xx = X 2 = A, r/ i = f?02 = ^, and Coi = C02 = C- V has 
a minimum at the same expectation value for $! and $2 and the Majoron coupling to j'J 
vanishes. 

Since the couplings to fermions are different for Hi and H 2 , the above exchange symmetry 
is no longer exact once fermion-loop corrections are included. If we adopt the "effective- 
potential" method for evaluating the quantum corrections to V, they are typically of order 
7 4 (7 is the Yukawa coupling of Hi to the neutral fermions). We assume these corrections 
are fine-tuned so that they will not have a significant contribution to the coupling. 

We now find out explicitly the particle spectrum in the Higgs sector. For the charged 
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scalars h a , the mass matrix is given by 

/ 2(C + ^K "(C + 0¥L -({ + Z)v v L 

-(( + ()¥l (C + eK 2 + Ci2^i -Ci2«i 



(9) 



The mass eigenstates and their masses are readily evaluated: 
h G 



ril sin 3 . , , . 9 

cosp/i H ^(^1 + ^2); m —0, 

v2 



v/2 



(h 1 -h 2 ); m 2 = (C + eK + 2 Ci2^, 
cos/9 



-sin^o + ^^i + M; ™ 2 = (C + OK 2 + 2 ^i), 



where the angle 3 is defined by 



tan/5 = 



\/2v L 



v 



(10) 

(11) 
(12) 

(13) 



The zero eigenstates h G correspond to the NG bosons associated with the SU(2) breaking 
and they are absorbed into the W 1 *" gauge bosons. For the real part of the neutral sector, 
the mass matrix is 

4A ug 2rjv v L 2rjv v L 
2rjv v L A\v 2 l + fug 2r] 12 v 2 L -fug | , 
2r]v v L 2r] 12 v 2 L - fug 4Au| + fug 

and the mass eigenstates and masses are 



(14) 



01 



(0i - 2 ); m 2 = (4A - 2 Vl2 )v 2 L + 2fug, 



V2 

cos a 0o + 



sin a . . o o 4 

—j=-((p 1 + (p2); m = 4A u + tana -j=r}v v L , 

cos a: , , , s 2 2 . J 



0- = -sina^oH — (01 + 02); m = (4A + 2r] U )v L - tan a -j=riv v L , 

V2 



where a is defined by 

cot a — tana 



V2 

V2 
V 



(15) 
(16) 
(17) 



A 


(A + f) 


Hi 


L v L 







7T 7T 
--<«<-. 

4 4 



(18) 



Similarly for the imaginary part, the mass matrix is 

/ 4u£ -2v v L -2v v L 
f -2v v L ug v 2 
\ -2v v L Vq ug 



(19) 
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The mass eigenstates are 



COS U 

^ A = -an/3 (po + —^-( ( p 1 + ( p 2 ); m 2 = 2£(v 2 + 2v 2 L ), (20) 
= -j=((pi - (p 2 ); m 2 = (21) 

Sill u 

(p G = cos(3<p + —^-(^ + (£2)] m 2 = 0. (22) 

The combination tpo is absorbed into the Z° gauge boson and (pL is the Majoron. 

The motion of the Majoron, especially its interaction with the other fields, is now solely 
described by the conservation, Eq. fl5|). We obtain the effective interaction Lagrangian 
£ e ff by requiring that the resulting Euler-Lagrange equation with respect to coincides 
with Eq. (|). Up to cubic terms, we get 



C 



cff 



. pt® + -^Z^ L -gv L (W+hl + W-h L 
vl / L2 cos9 w v p f 

+7: (d^if L ) (cos a 0_ + sin a 0+) + cos /3 (<9 m </?a) 4>l 



+ 



2 

2 



cos/3 - ^(^/i L ) + (^4)^ - 4(^M} 



(23) 



Reflecting on the NG-boson nature of <pl, we have written its interaction in the derivative- 
coupling form. The first term in the bracket includes the coupling to the neutrinos and 
induces neutrino decay. The second term is important since it describes the coupling to Z°. 
Remember that the mass of the accompanying scalar particle, 0l, is (4A — 2r]i2)v\ + 2£t>Q. 



The values for v and vl are constrained by yt>o + 2v\ = 174 GeV, but otherwise they are 
free. The expression for the 4>l mass can naturally give a bigger value than the Z° mass. 
Hence the decay Z° — ► (Pl4>l can be forbidden kinematically or else this model would be 
ruled out by the present LEP experiments. Note also that in the DM model all the Higgs 
bosons belong to the doublet representation and vl and Vq are free from the p-parameter 
constraint. 



8 



3. Neutrino decay 

We now look into the neutrino decay 

Ur-^u^cfL or V e Lf L . (24) 

They are induced by the first term in the bracket in Eq. (^). The existence of these flavor- 
changing processes can be seen in two steps. First, the conservation of the neutrino portion 
of j^(f) is violated in the symmetry-broken phase by the appearance together of both the 
Dirac mass, 

mBab = v L T\ b , (25) 

and the Majorana mass M^. This explains why the Majoron coupling to a charged lepton is 
suppressed even though j^(f) has charged- lepton components: their masses are necessarily 
of the Dirac type, hence that part of is automatically conserved at tree level. Second, 
although j^(f) is diagonal with respect to z/ a , the mass diagonalization procedure involves 
both v a (L = 1) and N a (L = 0), and this generates nondiagonal flavor- changing vertices. 
This is in contrast to the TM model, where no singlet neutrino is involved and the current 
is flavor-diagonal even after the diagonalization of the Majorana mass matrix for v. 

Let us see this second point in detail for the seesaw mass matrix and obtain the neutrino 
decay vertex explicitly. In the symmetry-broken phase the physical fields are the mass 
eigenstates. We write the neutrino fields in j^(f) in terms of these. It is done by the 
following replacement: 

UHU)-(££)U)- 

where we have suppressed the generation index; Ui are 3x3 matrices. The unitary matrix 
U is obtained by 



Although U is unitary, U\ is not necessarily so. Thus the light neutrino components in j^(f) 
after the replacement, u^UiUia^u, are not diagonal in general and allow the flavor-changing 
decay. For a typical seesaw mass matrix (mrj <C M), U2 and U3 are of order m^/M and the 
deviation of U\ from a unitary matrix is small, of order (mo/M) 2 . The matrices £/, have 
been obtained order by order in tud/M ||. f/i is given explicitly up to second order as 

Ui = (l - ^d^^E) V, (28) 

where the unitary matrix V is defined by 

V T ( _m D^ m 5) ^ = diag(m^, m^m^), (29) 

and in solving [/j we have assumed without loss of generality that M a b is already diagonal 
with positive eigenvalues. The eigenvalue m u has the typical seesaw size, mf^/M. The 
neutrino decay vertex is now given by 

^T^EVU^**), J^V+m^mSV (30) 

and the flavor- changing mixing is of order m v /M. 

The decay width of z/ T is readily evaluated. For m Vr ^> , m Ue , we obtain 

f V- 1 ip |2 m ^ 7 2 sin 2 g (m„\ 3 
af^e 64vr 64vr \ M J 

where we parametrize \R\ as (m Vl sin 9/M) with a mixing angle # and {m UT /vL) 2 as ( , y 2 m UT /M) 
with the Yukawa coupling 7 for the r neutrino. 

In the SM model this decay is also similarly described but the power of itl Vt /M in the 
final result ([H|) is 4 instead of 3. This is because the corresponding scalar expectation value 
vl in the SM model is related to the mass of the gauge-singlet neutrinos. 

Let us compare this width with the cosmological bound. The thermal history of neutrinos 

in the expansion of the Universe is well studiedF] |L8|. Their relic density parameter f2„ p is 

1 The key process to get the relic energy density of a neutrino species in Ref. |Q is vv «-> / / '. In the 
following, we assume that extra processes, intrinsic in a Majoron model such as v v «-> ip^ <Pl, will not change 
the results drastically. 
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given by 



for each light neutrino, i.e. that which is relativistic at its decoupling temperature, or by 



for each heavy neutrino, i.e. that which is nonrelativistic at its decoupling temperature. 
[The notation is the same as in Ref. |Tj|. Since we are interested in a narrow mass range, 



around 10 MeV, we have neglected a logarithmic dependence on m v in Eq. (^).j Thus a 
neutrino species that is heavier than 100 eV or lighter than a few GeV must decay. For 
the lifetime constraint we use here that which comes from a consideration on "structure 
formation" fflSf] , which gives a much stronger bound than the condition of not overdosing 



the Universe by the relic density of the decay products f20 |. This requires the relic density 
parameter of the decay product to be smaller than that of the radiation, i. e. the photon and 
the light neutrinos, 

n^h 2 ~4x 10" 5 . (34) 

The relic density of the relativistic decay product decreases faster than that of the nonrela- 
tivistic matter and is approximately given by Q u9 Rd, where Rn is a scale parameter related 
to the lifetime T^ 1 by 

! 4.9 x 10 9 ^ 

r ~ (o^)V2 year " (35) 



Thus we get 



< 3 x 10 2 ( -^-] sec for light neutrino, (36) 
VIMeV/ 

/ 777, \ 4 

r _1 < 4 x 10~ 2 ( — ) sec for heavy neutrino. (37) 

VIMeV/ v 1 

We demonstrate in Fig. 1 that it is indeed possible for v T to be as heavy as 10 MeV for a 
choice of parameters 7 2 = 10~ 2 and sin 2 6 = 10 -2 as an example. We depict the allowed area 
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in the (m^—M) plane. The lines AB and BC come from the cosmological bound for light 
and heavy neutrino, respectively. They are obtained by Eq. (|36| ) or ([37]) with (|31|). We have 
also taken into account a constraint coming from supernova cooling by Majoron emission ||, 
which is shown by the line CD. The region allowed by the cosmological constraint mostly 
corresponds to that of "Majoron trapping" and it gives H 

/ m v \ /lGeV\ 2 o . , 

Note that the values for m Ur and M implicitly determine the value of Vl by the relation 
jvl ~ \/ Mm VT . There is of course also the laboratory upper bound of 35 MeV on itl Vt . 

The values of for the allowed region in Fig. 1 are relatively low compared with vq ~ 
170 GeV; they are typically a factor of 5 or so less. This predicts a relatively light Higgs 
boson, the </>_ of Eq. (17). Since this Higgs boson is orthogonal to <pL which couples the 
Majoron tp^ to the Z° boson, there is no conflict with the experimental data for the Z° width. 
It can, however, be a rare decay product of Z° in 

Z ^ </>_//. (39) 

The branching fraction of this process is suppressed by 

(— cos/5 sin a + sin/5 cos a) 2 ~ 1CT 2 (40) 

compared with that of the single Higgs boson of the minimal standard model and thus 0_ 
can still have a mass below the latter's experimental lower bound of about 60 GeV. Through 
its (po component, </>_ decays into visible channels such as charged fermion pairs, and may 
thus be observed in future Higgs-boson search experiments. We have drawn the line AE that 
corresponds to Vl = 10 GeV in Fig. 1; 0_ with this mass has roughly the same branching 
fraction as the standard Higgs boson of 60 GeV in the process (^) plf . 

Finally we consider the effect of nucleosynthesis. The number of relativistic degrees of 
freedom at the temperature of about 1 MeV strongly affects the abundance of light elements 
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I8| . Using the best current nuclear physics data and astronomical observations, the number 
in terms of corresponding light neutrino species, N„, is restricted to be less than 3.3 [ fTU| . 
The Majoron contribution to N v has been considered extensively in Ref. ||. The key process 
for the Majoron to be in thermal contact with the neutrinos is vv <-> ipi^L- Its decoupling 
temperature Td is estimated by comparing the inverse mean free path of the process with 
the expansion rate of the universe, 

Td ~ 10V^(T D )^|— . (41) 
v mpi 

The decoupling temperature from the r neutrino is the lowest among the three neutrinos and 
it is much lower than 1 MeV for the parameters that are allowed in Fig. 1. The Majoron keeps 
thermal equilibrium as long as the r neutrino does. Thus its decoupling temperature cannot 
be high enough to suppress its contribution to N v [[J, which is given by (8/7) (1/2) ~ 0.6. 

To keep N v within the above-mentioned bound, the r neutrino needs to be nonrelativistic 
at its decoupling temperature and its energy density must be small compared to that of the 
radiation. Since the enery density of massive matter after decoupling decreases at a slower 
rate than that of the massless degrees of freedom, it may significantly contribute to the 
total energy density again afterwards. If this matter domination takes place at the time 



of nucleosynthesis, it would affect the primordial abundance of light elements p2| . The 
cosmological constraint on the r neutrino lifetime we have used does not allow its energy 
density or that of its decay products to exceed that of the light degrees of freedom. But it 
may still contribute partially to N v . Based on a detailed study of the r neutrino lifetime 
and mass constraints from nucleosynthesis, we thus further require its lifetime to be shorter 



than 1 second |22|, |23[. We show the corresponding boundary by a dashed line in Fig. 1, the 
region to the right of which remains allowed. Under the conditions that r neutrino is strongly 
nonrelativistic at the decoupling temperature and it decays fast enough before primordial 
nucleosynthesis, N v is given by the sum of the contributions from the Majoron and the two 
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lighter neutrinos, i.e. 2.6. 

In summary, we have presented a new Majoron model. It allows the r neutrino to be in 
the 10-MeV mass range. 
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Figure caption 

Fig. 1 Region of the (m VT -M) plane allowed as the result of various constraints for 
sin 2 # = 1(T 2 . 
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